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Abstract 

We analyse the structure of long distance (LD) contributions to the CP-violating parameter e^, that 
generally affect both the absorptive (T12) and the dispersive (-M12) parts of the K° — K° mixing 
amplitude. We point out that, in a consistent framework, in addition to LD contributions to ImF^, 
estimated recently by two of us, also LD contributions to ImMi2 have to be taken into account. 
Estimating the latter contributions the impact of LD effects on ex is significantly reduced (from 
—6.0% to —3.6%). The overall effect of LD corrections and of the superweak phase being different 
from 45° is summarised by the multiplicative factor n e = 0.94 ± 0.02. 

1 Introduction 

Some of the most important tests of the Standard Model (SM) are offered by CP- violating observables, 
that in this model are supposed to originate from a single CP-odd phase in the CKM matrix jl . In 
particular, the crucial test is the hierarchy of CP-violating effects in F>d, B 8 and K systems predicted 
by this model. Indeed the most prominent CP-violating observables in these three systems, S^,k S j >$W 
and €k, predicted by the SM, differ by orders of magnitude 

S^ Ks &2/3, « 4 x 1(T 2 , |e K |«2xlO- 3 . (1) 

Extensive analyses of the Unitarity Triangle have shown a spectacular consistency of the data for 
S-ipK s an d e K, within the parametric and theoretical uncertainties in ex, that until recently were rather 
sizable. The size of measured by CDF |2J and D0 [3| appears to be by one order of magnitude 
larger than predicted by the SM, but the large experimental errors preclude any definitive conclusions. 

Recently the consistency of the measured values for S^k s an d £k within the SM has been challenged 
in [ij[5j due to two facts: 

• The improved value of the relevant hadronic parameter Bk from unquenched lattice QCD that 
enters the evaluation of ejf. This parameter is now not only known with an accuracy of 4% [6j[7] 
but turns out to be significantly lower than previously found in lattice calculations, suppressing 
by 10% the previous estimates of ex- 

• A more careful look at e^, that identified an additional suppression of \ck\, summarised by a 
multiplicative factor n e — 0.92 ± 0.02 15] to the previously adopted formula for ck- 
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In view of these two suppressions, as demonstrated in |5], the size of CP violation measured in 
Bd — > ipKs might be insufficient to describe ex within the SM. Clarifying this new tension is important 
as the S^Ks ~ e K correlation in the SM is presently the most important direct relation between CP 
violation in the Bd and K systems that can be tested experimentally. 

The correction calculated in [5] originates from two factors: i) the difference of the superweak 
phase (f> e from 45° , and ii) the long-distance contribution to €k arising from the imaginary part of the 
absorptive amplitude of the K° — K° mixing, T^. The latter effect has been estimated with the help 
of the AI = 1/2 dominance in K — > tttt decays and the experimental value for e'/e. 

In the present paper we point out that at the same level of accuracy other effects should be con- 
sidered, in particular the long distance contributions to the imaginary part of the dispersive amplitude 
M 12 . While this topic has been the subject of intensive discussions in the mid 1980's, it is important 
to have a fresh look at this issue in view of the decrease of the error in Bk and of the theoretical 
advances during the last twenty five years. 

Our paper is organized as follows. In Section [2] we present general formulae from which the 
different contributions to ejf can be clearly identified. In Section [3] we discuss €k using the Operator 
Product Expansion (OPE). This allows us to identify the most important, still missing, long-distance 
contributions to ImM 12 . In Section [4] we estimate the size of these contributions in the framework of 
Chiral Perturbation Theory (CHPT), and briefly compare our findings with previous literature. We 
conclude in Section [5] 



2 Notation and general formulae 

Indirect CP violation originates in the weak phase difference between the (off-diagonal elements of 
the) Hermitian matrices M and Y which control the time evolution of a neutral meson system. For 
the K° — K° system one has 

±(\K°(t)) \ _/ HZ |*0(t)) N 

l dt[ |tf°(t)> -{ M l 2 { \K°(t)) J ' ™ 



Defining the eigenvectors 



]Ks{L)) = TWtW) [(1 + m ° } T{1 ~ ' (3) 



the following phase-convention-independent relation holds: 



Re(e) _ Im(r 12 M* 2 ) 



HO Im 



12 



M 



12 



(4) 



This represents indeed the indirect CP-violating parameter measured from the semileptonic charge 
asymmetries [8] or the Bell-Steinberger relation [9]. The experimental smallness of Re(e) makes the 
expansion to first non-trivial order in the weak phases an excellent approximation. At this level of 
accuracy we can identify Re(e) with the real part of the complex quantity 6k , defined in terms of the 
K — > 2ir amplitudes, 

2??+-+77oo A(K l ^tt 1 ^) 
(K = 3 ' ^ = A(K s ^n^) ■ (5) 

The two parameters are indeed related by €k = £ + where £ is the weak phase of the K° — > (2tt)j = q 
amplitude, namely 

C=^, A = A(K°^(2n) I=0 ). (6) 

Expanding to first non-trivial order in the weak phases we have 

ArriK = rriL - ms = 2Re(M i2 ) , 

Ar = Y S -Y L = -2Re(r 12 ) . (7) 
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Introducing also the so-called superweak phase, (f> e = arctan (2Amx/Ar), the expression for Re(e) 
becomes 

[ ImAf 12 imr 12 

Re(6K) = Re(e) = cos<p e sm 



2ReMi2 2Reri2 



(8) 



A further simplification arises by the observation that the |(27r)/ = o) final state largely saturates the 
neutral kaon decay widths. Since 



r 



21 



r 



the |(27r)/ = o) dominance in the sum over final states implies 



imr 



12 



Refis 



J Re^ 



-2£ • 



Expressing ReAfi2 in terms of Am^ and using Eq. ( 10 1 we arrive at 



Re(e/f) = cos(/> e sin< 



ImA/12 



which is consistent with 



sin< 



ImAfi 2 



(9) 



(10) 



(11) 



(12) 



The equation above allows us to calculate 6k by taking (f> e and ArriK from experiment and calculating 
ImM 12 and £ in a given model, in particular the SM. In Ref. [H] only short distance contributions to 
ImMi 2 , represented by the well known box diagrams, have been included, while £ has been calculated 
by relating it to the ratio e'/e and taking the latter from experiment. As we will discuss in the following, 
this approach is not fully consistent: in this way ImTi2 and ImA/12 are evaluated at a different order 
in the OPE. In particular, long distance contributions to ImA/12, which are of the same order of ImTi2 



(the latter giving rise to the £ term in Eq. (12)), are missing 



3 Decomposition of Re(e^) using the OPE 

As shown in Eq. ^ the evaluation of ex requires the knowledge of the weak phases of both Mi 2 and 
Ti2. In this respect, we should emphasize that ImMi2 and ImTi2 are both generated at 0(G 2 F ). Since 
ReA/12 and ReFi2 are very similar in size (<j> e « 43.5°), we should consistently evaluate ImMi2 and 
ImTi2 at the same order in the OPE. 

The relevant effective Hamiltonians are 'Has=2 (contributing to ImA/ 12 only) and Has=i (con- 
tributing to both ImM 12 and Imr 12 ). The leading term in the OPE is the short-distance contribution 
to ImMi2, 

' * (<^°i?4 6 L 2 i^ })* (13) 



2 m k 



-Im 1 



where 



H 



(6) 

AS=2 



= (sd) v _ A (sd) V - A , 



(14) 



is the dimension-six AS = 2 effective Hamiltonian. The operator does not mix with other 

operators and the imaginary part of its Wilson coefficient is dominated by terms proportional to the 
top-quark Yukawa coupling^ At this order in the OPE one is neglecting terms generated by two 
insertions of AS — 1 operators (see Fig. [I]) which cannot be absorbed into the coefficient of . 
For consistency, this implies one should set ImTi2 to zero, since ImFi2 is the absorptive part of the 



1 The explicit expression of the coefficient function Fq, depending on quark masses and CKM elements, can be found 



in 10 . 
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Figure 1: Contractions of the leading |A5| = 1 four-quark effective operators contributing to M\i at 



diagrams in Fig.[T] In other words, the leading order result is obtained with the following substitutions 
in Eq. |TT]): 

ImMia -> IxnM$ = ImMf 2 D and £ ^ . (15) 
Going one step forward requires taking into account: 

1. non-local contributions to both ImM 12 and Imr 12 generated by the 0(Gp) dimension-six AS = 1 
operators, 

2. local contributions to ImMi 2 generated by dimension-eight AS* = 2 operators of (D(G 2 F ). 

The structure of the subleading terms in ImMi2 is very similar to the 0(G 2 F ) long-distance contribu- 
tions to K — >• ttvv, discussed in Ref. [IT]. The relevant effective Hamiltonian changes substantially 
if we choose a renormalization scale above or below the charm mass. Keeping the charm as explicit 
degree of freedom, dimension-eight operators are safely negligible and the key quantity to evaluate is 



d 4 x(K°\T 



1AS| = 



1AS| = 



i(0) \K°) 



(16) 



where the superscript in %as=i denotes that the we have two dynamical up-type quarks. The ab- 
sorptive part of 71 2 contributes to Ti 2 , while the dispersive part contributes to M\%. In the latter case 
the leading term in the expansion in local operators should be subtracted, being already i nclu ded in 



lmMi2 • I n principle, extracting the subleading contribution to ImMi2 directly from Eq. (|16[) is the 



best strategy: the result would be automatically scale independent. However, in practice this is far 
from being trivial also on the lattice, given the disconnected diagrams in Fig. [T] 

Following a purely analytical approach, we can integrate out the charm and renormalize %as=i 
below the charm mass. This allows to identify £ with the weak phase of the Aq amplitude, that, as 
mentioned, has already been estimated in Ref. [5] (see also [l2]). On the other hand, ImMi 2 assumes 
the form 



ImM 



12 



where ImM™ 11 ' 



-local 



ImMf 2 D 
r(8) 



ImAf™ 11 



-local 



(8) 



+ ImMj 2 



(17) 



and ImM^ are not separately scale independent. The structure of the dimension- 
eight operators o btai ned integrating out the charm, and an estimate of their impact on €k , has been 



presented in Ref. 13 . According to this estimate, ImM 12 is less than 1% of the leading term 



(81 

The smallness of ImM 12 can be understood by the following dimensional argument. First, it should 
be noted that the CKM suppression of the dimension-eight operators is (V* s V c d) 2 , namely the same 
CKM factor of the genuine charm contribution in r H^fg_ 2 - Second, even if we are not able to precisely 
evaluate the hadronic matrix elements of the dimension-eight operators, we expect 



(K°\Q[ 8) \K°) = 0(1) xm 2 K x (K°\QW\K 



(18) 



According to this scaling, the contribution of ImAfj^ is an 0{m 2 K /m 2 sw 15%) correction of the 
charm contribution (charm-charm loops) to ImM 12 \ which itself is an 0(15%) correction of the total 
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dimension-six contribution. We are thus left with an overall 0(2%) naive suppression of ImM^ with 



respect to ImAfj^. According to the explicit evaluation in Ref. 
even smaller. 

The only potentially large long-distance contribution to ImM 12 is the contribution of the non-local 
terms enhanced by the AI = 1/2 rule. For this purpose, we observe that if we had a single weak 
operator in Has=i, this would generate the same weak phase to both ImMf^ and IruT^. As we 
discuss in more detail in the next section, this is what happens to lowest order in CHPT, where the 
A J =1/2 part Has=i has only one operator, with effective coupling G$. Decomposing lmM^P as a 
leading term proportional to G§, and a subleading term with different effective coupling 



13 



the actual numerical impact is 



ImMff = ImMff I 



ImMh D \ 

LA Inon— Gi 



(19) 



we can write 



ImM^ 2 = ReMia 



lG s V1 I OS Re[(G^) 2 ] ' 

and identify the weak phase of Gg with £. As a result, 



ImM^I^ « ReMf^l^ x (-20 « -£ x (Ato^. 



(20) 



(21) 



This allow us to re- write Eq. (11) as follows 



Re(eif) = cos </> c sin < 



ImM- 



(6) 
12 



ArriK 



3lmMi 2 



(22) 



where <5i m M 12 encodes the subleading terms in ImM ] ^ £> | non _ G 2 (including also ImM^'). Note that, in 
the limit where the contribution of Gg saturates Arajc, the contribution of £ would be absent. This is 
exactly what we should expect, since in this limit M\% and Ti2 would have the same weak phase but 
for the short-distance contribution to ImMi2. 



4 Estimate of long-distance effects in CHPT 

A convenient framework for estimating the long-distance contribution to M\% is provided by Chiral 
Perturbation Theory (CHPT). In this framework 7r, K and r\ fields are identified with the would-be 
Goldstone bosons arising from the SU(3)l x SU(3)r — > SU(3)l + r symmetry breaking of the QCD 
action in the limit of vanishing light quark masses. Low-energy amplitudes involving these mesons, 
expanded in powers of their masses and momenta, are evaluated by means of an effective Lagrangian 
written in terms of the pseudo-Goldstone boson fields. 

The lowest-order effective Lagrangian describing non-leptonic AS — 1 decays has only two opera- 
tors, transforming as (8^, 1r) and (27 l, 1r) under the SU(3)l X SU(3)r chiral group. Moreover, only 
the (8^,1/j) operator has a phenomenologically large coefficient, being responsible for the enhance- 
ment of AI —1/2 amplitudes. As a result, the only term in the effective Lagrangian relevant to our 
calculation is 

(23) 



where, as usual, we define 



U = exp(iV2<S>/F) , 



$ = 



V2 



n 



7T 



7T+ 
/ , JL 

K° 



K+ 

K° 

_ *L 

V6 



(24) 
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Figure 2: Tree- level and one-loop diagrams contributing to K°-K° mixing in CHPT. 



and F can be identified with the pion decay constant (F sa 92MeV). The effective coupling G$ can 
be determined by K — > 2ir amplitudes. Neglecting the (27l, Ijj) operator and evaluating the K — > 2ir 
amplitudes at tree level leads to 

A Q = A(K° -> (2vr) /=0 ) = V2FG s (m 2 K - mj) , (25) 

which implies \G 8 \ s» 9 x 1CT 6 (GcV)~ 2 . As far as the weak phase of G 8 is concerned, at this level of 
accuracy we have Im(G8)/Re(Gg) = £. 

In principle £[as|=i cou ld contribute to M\% already at 0(p 2 ), via the tree-level diagram in Fig. jij 
(left). However, considering the 0(p 2 ) relation among 7r°, 77 and kaon masses (i.e. the Gell-Mann- 
Okubo mass formula), this contribution vanishes |14| . As a result, the first non- vanishing contribution 

to M12 generated by £^2s\=i ar i ses oru y at C(p 4 )- 

At 0(p 4 ) we should evaluate loop amplitudes with two insertions of £[as|=i an< ^ tree-level diagrams 
with the insertion of appropriate 0(p 4 ) counterterms. Among all these 0(p i ) contributions, the only 
model-independent, and presumably dominant, contribution to M12 is the non-analytic one generated 
by the pion-loop amplitude in Fig. [2] (right), 

3 ^(g; 



T-, 



12 



y/1 - 4r 2 f log 



16tt 2 



J2\1 



K 



1 + y/1 - 4r 2 
l-^/T^rf 



log 



(26) 



with r 2 . = rn 2 ./m 2 K and where we have absorbed all finite (mass- independent) terms in the definition 
of the renormalization scale \x. This is the only contribution which has an absorptive part. As a 
consequence, its weak phase can be unambiguously related to the weak phase of the K° —> (2tt)j =0 
amplitude to all orders in the chiral expansion. In addition, it is the only contribution that survives 
in the limit of SU(2)l x SU(2)n CHPT, which is known to represent a good approximation of the 
full 0(p 4 ) amplitude in several X-decay observables where contributions from counterterms are fully 
under control (see e.g. Ref. [15]). 

A CHPT calculation of M 12 complete to 0(p 4 ) would require consideration of loops involving kaons 
and 77's, as well as 0(p A ) local counterterms. However, all these additional pieces are not associated 
with any physical cut. As such, they can effectively be treated as a local term whose overall weak 
phase cannot be related to the phase of the K° — > (2tt)j =0 amplitude^] On account of the above 
considerations]^] we refrain from a full 0(p 4 ) CHPT calculation, and we focus on the pion-loop non- 



analytic contribution only. Using the relation T^"' — 2ttikM{! 1 (/z), the result in Eq. (261 implies 



m[y\^) = - 



* \2 



647r 2 m^- 



log 



' K 



o 



K 



(27) 



The absorptive part in Eq. (26 1 is nothing but the leading |(27r)/ = o) contribution to Pi2, which gives 
rise to the relation (10 1. The dispersive part is the dominant contribution to M 12 in the leading-log 

2 For a recent, elucidating discussion about the role of kaon loops in CHPT, see 16 . 

3 The authors warmly acknowledge Jean-Marc Gerard for triggering a discussion on this point. 
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approximation. The close link of these two terms is a further confirmation that we cannot neglect the 
long-distance contribution to ImMi2 if we want to keep track of all the 0(£) terms in ck ■ 

Using the result in Eq. (27) we can estimate the contribution to ImMi2 proportional to G$, which 
enters in the phenomenological formula for Rc(eA-) in Eq. (22 1. Setting p = 800 MeV and varying it 
in the interval 0.6 ~ 1 GeV leads to 



Am L K D \ 



A) 



cxp 



2ReikT 



(7T7T) 

12 



Aj 



cxp 



0.4 ±0.2 



(28) 



Note that the result has a well-defined sign since G$ (or ^4o) appears squared in ■ Using this 

result in Eq. ( 22 ) we find a suppression of the £ term relative to the estimate in [5] , where only the 
LD contribution to ImT]^ has been taken into account. 

Since our estimate of Amj^| G 2 is not the result of a complete calculation at fixed order in the 
chiral expansion, it is worthwhile to cross-check it using a different argument. For this purpose, we 
note that the only relevant contribution to M12, beside the two-pion intermediate state, is expected 
to arise from the tree-level 77' exchange (Fig. [2] left) [17]. We can thus decompose Myi as follows: 



Mi 



M 



SD 



M LD \ 

1V1 12 |W7T 



M LD \ , 

1V1 \2 W 



(29) 



According to this decomposition it is clear that, as far as long-distance contributions are concerned, 
we can trade the evaluation of M^f\ V7r for that of Mfifl^i. An estimate of the 77' contribution to M 12 
goes beyond pure CHPT, where it can be considered as a free parameter (the leading contribution 
to the 0{p A ) local terms). However, its impact can be estimated in the large N c limit, extending 
the underlying symmetry from SU(3)l x SU(3)r to U(3)l x U(3)r. Within this framework the 
operator basis must be extended and we cannot directly relate the phase of the rf exchange amplitude 
to the phase of G 8 . According to the recent analysis in Rcf. 17 , the 77' amplitude gives a negative 
contribution to Amu- 



2ReM[ 2 D \, r 



= Ai 



-0.3A? 



cxp 



(30) 



Most important for our analysis, this contribution is found to be induced at the quark level by the 
operator (sd)v-A x (uu)v-A only 17 . This implies that the 77' exchange has a vanishing weak phase 



in the standard CKM phase convention: 



ImMf/V = 



(31) 



Using this result in Eq. (11 1, and using the relation (21 ) for the tttt contribution, we get 



R-e(e^) = cos</> £ sin< 



ImM- 



(6) 
12 



cxp 
l K 



Am S pP + Am%P\ n > 
A^fP 



(32) 



where the Gs term (i.e. the tttt contribution), is manifestly absent. Denoting as p the coefficient of the 



AmfP 



(0.7±0.1)Am£ r 
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£ term in Eq. (32 1, and combining Eq. (30) with the NLO short-distance estimate of ReMi2, namely 



we get p = 0.4 ± 0.1. This result is well consistent with the value 



0.6 ± 0.2 obtained from Eq. ( 22 ) with the direct evaluation of the tttt contribution in Eq. ( 28 1 



We rate the direct evaluation of the tttt loop as the most reliable estimate of p. As a consequence, 
our final phenomenological expression for €k is 



€k = sin ( 



ImM 



(6) 
12 



ArriK 



with 



p = 0.6 ±0.3 



(33) 



where we have conservatively increased by 50% the error in Eq. ( 28 ) to take into account the sub-leading 
contributions of ImM ] Z 2 £> | non _ G 2. For p — 1 our result reduces to the one in 5 . The contribution 
calculated in this paper, resulting in p < 1, completes the estimate of the terms of 0(£) in ex ■ 
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Following the notation of Ref. [5], we summarise the corrections to due to LD effects and 
7^ 45°, via the introduction of the phenomenological factor n e , defined by 



AniK 



(34) 



According to our result in Eq. ( [28] ), and taking into account the estimate of £ obtained in jSJ, namely 
£ = -(6.0 ± 1.5) x 10~ 2 x y/2\e K \, the new numerical value of n e is 

x [l + p-f =f — =0.94 ±0.02. (35) 



1/V2 V V2| 



This should be compared with 0.92 ± 0.02 in (5 and 0.92 ± 0.01 in 19 , where only the long-distance 
contributions to ImFi2 (not those to ImMi2) have been included. 

4.1 Comparison with previous literature 

As anticipated in the introduction, the relative role of short- and long-distance contributions to ex 



has been widely discussed in the literature in the mid 1980's 20 -28 . It is therefore useful to compare 
our findings to those in these earlier works. 

First of all, we agree on the main conclusion of all these papers, namely that e K D /e^ is small as 
long as e' /eK is small. This is certainly correct, but it is not the point of our analysis: the issue we are 
addressing in this work is the size of the subleading (long-distance) contributions to ex, that vanish 
in the limit of vanishing e' . 

Second, we agree that single-particle intermediate states (tt , ry, rj') do not generate a significant 
long-distance contribution to ImAf 12 . The cancellation of 7r° and r\ contributions at the lowest order 



in the chiral expansion was noted first in 24 . The role of the rf was more debated [24 -27 . The 
issue was clarified in |28| , where it was shown that the full nonet contribution (ir°, 77, ij') vanishes in 
the large N c limit. This is consistent with our findings, which are based on the updated and detailed 
analysis of the 77' exchange amplitude in Ref. [17] . 

Having clarified that single-particle intermediate states do not generate a significant contribution 
to ImM^ 2 D , we are left with the two-pion intermediate state as the potentially leading contribution to 
ImAff 2 D . A naive estimate of this contribution at the partonic level seems to indicate that it is totally 
negligible; however, as we have shown, this is not the case because of the A J =1/2 enhancement of 
K — > 2-7T amplitudes. Our key observation is that, thanks to chiral symmetry and to the AI = 1/2 
dominance, the weak phase of this contribution can be related to £, and the problem is shifted to the 



evaluation of the two-pion contribution to Am^, as summarised in Eq. (21 1. The numerical impact 



of this contribution is then estimated in two ways: i) a direct computation of the tttt loop in the 



leading- log approximation, Eq. ( 28 ) , which provides a definite sign for this term; ii) the difference 
between the experimental value of Amx and the sum of its short-distance contribution and the other 
large long-distance contribution provided by the rf exchange, which allows us to perform the useful 
cross-check: 

Am L K D \ Gl « Am^ - [Am s K D + Am K D \ n >] . (36) 

We finally note that our estimate of the 0(£) corrections to ck is based on the dominance of the 
A J =1/2 amplitude in K — > 2tt decays. Given the experimental smallness of AI = 3/2 transitions, 
and the overall size of the effect we have evaluated (a few % correction to ex), this is certainly a very 
safe approximation. 

5 Conclusions 

In this paper we have presented a complete analysis of £k beyond the lowest order in the OPE. In 
particular, we have analysed the structure of long distance (LD) contributions that affect both the 
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absorptive (r 12 ) and dispersive (-M 12 ) parts of the K a — K° mixing amplitude. We have pointed out 
that, in a consistent framework, in addition to LD contributions to ImTi2, estimated recently in j5j, 
also LD contributions to ImMi2 have to be taken into account. Estimating the latter contributions in 
chiral perturbation theory, we found that they reduce by 40% the total impact of LD corrections on 

The overall multiplicative factor n e in e^-, summarising the effect of LD corrections and of the 
superweak phase being different from 45°, is increased to n e = 0.94 ± 0.02, to be compared with 
0.92 ± 0.02 obtained without LD contributions to ImMi2. 
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